Lattice-induced non-adiabatic frequency shifts in optical lattice clocks 
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We consider the frequency shift in optical lattice clocks which arises from the coupling of the 
electronic motion to the atomic motion within the lattice. For the simplest of 3-D lattice geometries 
this coupling is shown to only affect clocks based on blue-detuned lattices. We have estimated the 
size of this shift for the prospective strontium lattice clock operating at the 390 nm blue-detuned 
magic wavelength. The resulting fractional frequency shift is found to be on the order of 10 -18 and 
is largely overshadowed by the electric quadrupole shift. For lattice clocks based on more complex 
geometries or other atomic systems, this shift could potentially be a limiting factor in clock accuracy. 

PACS numbers: 37.10.Jk, 06.30.Ft 



In 2003, Katori et al. [l[ first proposed a frequency 
standard based on neutral atoms confined in a well engi- 
neered optical lattice, citing the potential to combine the 
benefits of the state-of-the-art clocks based on single-ions 
in a Paul trap (Doppler-free and recoil-free spectroscopy, 
long interrogation times) and ensembles of neutral atoms 
in free-fall (large statistical samples). Over the last few 
years these attributes have largely been realized, with 
optical lattice clocks being demonstrated at 1 x 10~ 16 
fractional uncertainty @] ■ On a theoretical basis it is an- 
ticipated that these clocks will further reach to the 10~ 18 
level 0, H[ ■ To achieve such an unprecedented accuracy 
it is imperative that all frequency shifts induced by the 
lattice laser field are well understood and evaluated at 
this level. A number of important theoretical works have 
been dedicated to determining or estimating frequency 
shifts arising from mechanisms such as hypcrpolarizabil- 
ity (second order in intensity) [HH, black-body radiation 
o| , collisions Q , off- resonant couplings to the probe laser 
8| , and magnetic dipole and electric quadrupole interac- 
tions with the lattice field [1, @ . 

All optical lattice clocks which have been demonstrated 
to this point operate at lattice wavelengths which are de- 
tuned to the red of a nearby atomic resonance. Proposed 
clocks based on blue-detuned lattices are presumed to be 
inherently less susceptible to lattice-induced frequency 
shifts [lfj. This assumption is based on the fact that 
atoms in a blue-detuned lattice are localized to regions 
of intensity minima rather than intensity maxima, thus 
suppressing higher order hyperpolarizability effects and 
scattering of lattice photons. Takamoto et al. [l(j dis- 
cussed the prospects of a Sr clock with a blue-detuned 
lattice and have experimentally identified an optimal 
magic wavelength at 390 nm. Suitable blue magic wave- 
lengths are also likely to exist in other alkaline-earth(- 
like) species (e.g., Mg, Ca, Yb). To exemplify the com- 
mon conviction that blue-detuned optical lattice clocks 
will ultimately provide superior accuracy to their red- 
detuned counterparts, we refer to recent Letter Q in 
which it is argued that clock uncertainty at the 10~ 18 
level will be within reach; we note that this conjecture 
is preceded by the caveat of operation at a blue-detuned 



magic wavelength. 

In their original proposal, Katori et al. suggested the 
possibility of using 3-D lattices with single atom occu- 
pancy (or less) at each site, thus effectively eliminating 
collisional effects. Recently a 3-D lattice clock has been 
realized with bosonic 88 Sr atoms [l4j|, marking another 
important step towards achieving the anticipated accu- 
racy of the lattice clocks. 

Here we consider the frequency shift in optical lattice 
clocks arising from the coupling of the electronic mo- 
tion to the atomic motion within the lattice; we refer 
to this as the non-adiabatic shift. For the simplest 3- 
D lattice geometry we find that this coupling only ren- 
ders a frequency shift in the blue-detuned lattice clocks. 
We estimate the size of the non-adiabatic shift for the 
prospective "blue" Sr lattice clock and find it to be on 
the order of 10 -18 fractionally. In this particular case, 
the non-adiabatic shift is largely overshadowed by the 
electric quadrupole shift, though for other geometries or 
for clocks based on other atomic species it could perhaps 
provide the leading source of clock uncertainty. 

We begin our analysis with the total Hamiltonian for 
an atomic system in the presence of some external laser 
configuration (e.g., an optical lattice), 

H(r,t,q)=T(r) + H at (q)+H h (r,t,q), 

where r denotes the atomic center of mass and q en- 
capsulates all electronic spatial coordinates relative to 
the center of mass as well as spin coordinates. Here 
T(r) = — (fi, 2 /2A/)V 2 represents the kinetic energy as- 
sociated with the center of mass motion (M being the 
mass of the atom), H at (~q) is the conventional field- 
free atomic Hamiltonian (including electronic kinetic en- 
ergy, Coulomb interactions, Breit interactions, etc.), and 
ilL(r,i, q) represents the interaction of the atom with 
the laser field. We treat the laser field classically, and 
as such the operator Hi,(r,t,q) is a simple function of r 
and t (i.e., it does not contain derivatives with respect to 
these coordinates) ; furthermore it is assumed to be peri- 
odic in time with period r, H^(r,t + T,q) = H^(r,t,q). 
It follows that the Hamiltonian H(r,t,q) is also periodic 
in time with period r. 
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Due to the periodicity of the Hamiltonian, we may 
treat this problem within the framework of the well- 
known Floquet formalism [ll[, in which the evolution 
of the system is governed by the eigenvalue equation 

H(r,t,q)iP(r,t,q) = £i>(r,t,q). (1) 

where H(r,t,q) = H(r,t,q) — ih-^ and ip{v,t,q) is peri- 
odic in time with period r. With the total wave func- 
tion given by ^(r,t,q) = ->jj(r,t,q)e~" l£t ^ h , the eigenvalue 
equation (TTJ) is entirely equivalent to the time-dependent 
Schrodinger equation. The eigenvalue £ corres pon ds to 
the time-averaged energy of the state \I>(r, [12| and, 
for an optical lattice clock, therefore represents the phys- 
ically significant energy level being probed. We will refer 
to H(r,t,q) as the steady state Hamiltonian and Eq. ([T]) 
as the steady state Schrodinger equation. 

We wish to find approximate solutions to Eq. ([T]). To 
this end, we begin by considering the Hamiltonian exclu- 
sive of the kinetic energy term T(r). In particular, we 
look for solutions to the eigenvalue equation 

n at ^(r,t,q)x n (r,t,q) = W„(r)x„(r, t, q), (2) 

where the atom-laser steady state Hamiltonian is 
"Hat,L(r-,£,(7) = T-L(r,t 7 q) — T(r) and we have introduced 
the quantum number n to label the eigcnsolutions. Effec- 
tively the atom-laser steady state Hamiltonian describes 
an atom and its interaction with the laser field when the 
atom is at a fixed location r within the laser field. Be- 
low we will consider solutions to Eq. @ in the context 
of steady state perturbation theory; for the time being 
we will assume the eigenvalues and eigenfunctions are 
known. We will, however, find it useful at this point to 
separate the eigenvalues into two terms, 

W n (r) = W n + U n (r), (3) 

where W n is the energy of the atom in the field-free limit 
(Hi,(r,t,q) —> 0) and U n (r) is then the corresponding 
level shift induced by the presence of the laser field. 

As the second piece of our approximation we consider 
the motion of the atom within the laser field as governed 
by the time-independent eigenvalue equation 

[T(r) + U n (r)] ip nv (r) = G nv (p nv (r), 

with the eigensolutions (for a specific atomic state n) 
being labeled by the quantum number v. We assume 
these eigensolutions to be known as well. 

We now graft our decoupled eigenfunctions together, 

tp„v(r,t,q) = Xn(r,t,q)ip nv (r). 

With the approximation S7 2 Xn( r ,t,q)(p nv (r) = 
Xn(r, t, q)V 2 <p m ,(r), it subsequently follows that 

H(r,t,q)ip nv (r,t 7 q) = (W n + G nv ) ip nv (r,t,q), 

i.e., the ip nv (r,t,q) represent approximate eigenfunctions 
of the steady state Hamiltonian, Eq. ([1]), with corre- 
sponding eigenvalues (W n + G nv ). Generally speaking, 



for atoms trapped in an optical lattice this approxima- 
tion is quite good and can be justified by the fact that 
the electron mass is much smaller than the nuclear mass. 
In physical terms, as an atom moves within a lattice site 
it experiences a varying lattice field and the electronic re- 
sponse to this varying field occurs on a time scale much 
faster than the atomic motion (the atomic motion effec- 
tively being adiabatic). We will refer to the correction 
which accounts for coupling of the electronic motion to 
atomic motion the non-adiabatic correction. We empha- 
size that the non-adiabatic effects discussed here are not 
associated with any type of dynamical change in the op- 
tical lattice field (e.g., ramping up or down the intensity). 

We aim to estimate the size of the non-adiabatic cor- 
rection. Analogous to the time-independent Schrodinger 
equation, the steady state Schrodin ger equation is deriv- 
able from a variational formalism [11], and as such we 
may estimate the true value of £ with an appropriate trial 
function, namely the approximate ( "adiabatic" ) eigen- 
functions ip nv (r,t, q). A reasonable estimate of the non- 
adiabatic correction, F nv , is then given by (we assume 
the adiabatic eigenfunctions to be normalized) 

F nv = - (W n + G nv ) , (4) 

where the somewhat cumbersome yet useful bra-ket no- 
tation with triple-angled brackets indicates integration 
over r, t, and q, with the i-integration over one period, 

(MUM)) = \j J + ^ J d 3 rdtdq 

XtP*(r,t,q)H(r,t,q)iP(r,t,q). 

We will also find it useful to designate the double- 
angled brackets ((. . . )) for integrations over t and q and 
single-angled brackets (...) for integrations over q alone. 
Eq. (Q| reduces to the form 

Fnv = ~&i J ^ [vL(r)<Pn V (r)((Xn\V 2 \Xn))(r) 

+2< w (r)V^(r)-(( Xn |V|xn»(r)], (5) 

where the notation ((x„|V 2 |x„))(r) and ((% n | V|x„))(r) 
emphasizes the fact that these are functions of r. 

Ultimately we are interested in the non-adiabatic shift 
of the clock frequency. The clock probes the energy in- 
tervals 

A£ v = AW + AG V + AF V , 

where A, here and below, is used to denote the differ- 
ence between clock levels, e.g., A£ v = £p.„ — £$, v for 
clock states 3 Pq and 1 5o, and we assume the usual case 
of v — > v transitions. Here AW is the unperturbed in- 
terval and the combination AG V + AF V gives the lattice- 
induced perturbation leading to a shift in the measured 
clock frequency. 

To estimate the size of the shift AF V , we return to the 
atom-laser steady state Schrodinger equation, Eq. ©. 



3 



Eigcnsolutions of Eq. ([2]) can be obtained by a systematic 
application of steady state (Floquet) perturbation theory 
[ill ] , treating the interaction Hl (r, t , q) as a perturbation. 
The zeroth order atom-laser steady state Hamiltonian, 

^itj, = ^at,L(r,t,g) - H L (r,t,q), has eigenvalues 
and cigenfunctions satisfying 



(6) 



We note the absence of r in this last expression. The 
zeroth order eigenvalues and eigenfunctions are given by 

W^=W n , X £\q) = fn(q), 

where W n and f n (q) are the eigenvalues and eigenfunc- 
tions of the time- independent Schrodingcr equation for 
the field- free atom, H at (q)f n (q) = W n f n (q). 

A complete set of eigensolutions to Eq. (JBJ) (spanning 
the space of all functions square integrable over q and 
periodic in time with period r) is provided by 



r v n p 



X { $(t,q) 



fn(q)e 



ipujt 



(7) 



where uj = 2tt/t and the additional quantum number p 
takes the value of any integer. The various values of p do 
not correspond to physically distinct states (i.e., the wave 

function Xnp (t,q)e~ lVVn v t / h is independent ofp), and we 
obtain a one-to-one correspondence with atomic states 
by setting p = above. However, in the application of 
steady state perturbation theory we require summations 
over a complete set of states, and this is furnished by 
Eqs. (0 with all integer p. We note that these are the 
familiar "dressed" atomic states (in which case p is simply 
associated with photon number). 

At this point it is useful to consider a specific form for 
the interaction Hjj(r,t,q). Initially we will limit our at- 
tention to the dominant electric dipole (El) interaction. 
The El interaction effectively "samples" the electric field 
at the location of the atom and has the form 

tf L (r,*,l) = -Re[E(r)e- iwi ] -T>(q), 

where D(g) is the atomic dipole operator and 
Re [E(r)e -1 "*] is the (classical) electric field with an as- 
sumed harmonic time dependence and complex ampli- 
tude E(r). 

The application of steady state perturbation theory 
proceeds in the same manner as conventional time- 
independent perturbation theory. The first order correc- 
tion to the eigenvalue, W^'(r) = ((Xn\Hi,\xn)), van- 
ishes when the integration over time is performed (and, 
within the El approximation, from parity considerations 
as well). After integrating over time, the second order 
correction reads 



M 2) (r) = j£ 



[E*(r) 



IDI 



[E(r) 



(8) 



where |n) = |/„) and the notation {. . . }® indicates that 
one should add to the term in braces a complementary 
term acquired by making the simultaneous substitutions 
E(r) — > E*(r) and oj -> —u. We will not consider 
higher order (hyperpolarizability) corrections, and so in 

line with Eq. ©, we set U n (r) = wl?(v). 

For scalar atomic states (e.g., clock states 1 S'o, 3 Pq) 
Eq. {5]) may be recoupled [13J and expressed as 



U n (v) = -a n I(r), (9) 

where I(r) = (c/8ir) |E(r)| 2 is the intensity at r and a n = 
(c/27r)a„ is the conventional (scalar) ac polorizability, 

© 

, (10) 



*»=4^ |(n ' |D|n>|2 { 



i 



with c being the speed of light. We note the in- 
dependence of the ac polarizability. 

After integration over time, the first order correction 
for the eigenfunction %„(r,t,g) is found to be 



ti ) (r,t,q) = -lY,{\ 



[E(r)>'|D|n)] t 



fw(q). 



For some operator S(r), the leading correction to 

((Xn\S\xn))(r) is given by ((Xn l^lxn ))(r). In par- 
ticular, we are interested in the cases S(r) = V 2 and 
S(r) = V. We now assume that the lattice is composed 
of a superposition of standing waves of linear polariza- 
tion, with all standing waves in phase (e.g., the "folded" 
3-D optical lattice in Ref. [l4[ ) ■ In such a lattice arrange- 
ment E(r) may be taken as real and we find the following 
expressions for scalar atomic states after recoupling, 



((xl 1) |V 2 |x«»(r)=-(^) 2 an/(r) ! 
((xl 1) |V|xl 1) »(r) = i ( T, i V/(r), 



(11) 



where the frequency-dependent atomic factor a n = 
(c/2Tr)a n is given by 



a„ = ^|MD|n)| 2 { J 

6 ^f \ (W n - Wn' 



(12) 



We now assume a specific lattice geometry; we choose 
the most straightforward 3-D geometry (case (I) of 
Ref. [|[) in which the lattice is composed of three stand- 
ing waves aligned in the x, y, and z directions having 
mutually orthogonal linear polarizations. It follows that 
the intensity is given by 



7(r) 



_ \l [4 cos 2 



IyCOS 2 



y 



TrCOS 2 



4 [I x sm 2 (^x) + IySin 2 (*y) + J z sin 2 (±z) 



0] 



(13) 

where have chosen the coordinate origin separately for 
the cases of red (top) and blue (bottom) lattices such 
that U n (r) (Eq. ©) has a minimum at the origin. 
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We now Taylor expand Eq. (|T3|) to quadratic terms 
in x, y, and z and evaluate the non-adiabatic correction 
(Eq. (0) using Eqs. (JTTJ) and assuming harmonic oscil- 
lator vibrational functions tp n v{r) with v = (v x ,v y ,v z )] 
the resulting expression reads 

op3/2=. 

F nv = mna n i + (-i)<+if£L_^ V ^ («, + 1) , 

(14) 

where I = (I x + I y + I z )/3, R = h 2 (uj/c) 2 /2M is the 
recoil energy, and t is equated to and 1 for red and 
blue lattices, respectively. 

To understand the influence of F nv on the clock fre- 
quency, it pays to consider an explicit form for the vi- 
brational energies G nv as well. Taichenachev et al. Q 
have shown that magnetic dipole (Ml) and electric 
quadrupole (E2) interactions can give non-negligible con- 
tributions to G nv , and so we include these effects in our 
consideration. For the present lattice arrangement, G nv 
is given by Q 

G nv = -I2a' n h + WR l«nl Yl ( Vi + i) ' ( 15 ) 

i 

where again we have employed the harmonic oscillator 
approximation. Here a n now includes Ml and E2 ef- 
fects in addition to the dominant El effects (a n = 
a^ 1 — a^ 11 — a^ 2 ) and a' n is the small corrective dif- 
ference (a' n = a„ + a;^ 2 ); this slightly modified a„ is 
appropriate in Eq. (jT4j) as well. We note the simlari- 
ties between G nv and F nv : both expressions consist of 
two terms, with the first term being independent of vi- 
brational state and proportional to Iq, while the second 
term is proportional to J^i VF ( v i + h)- 

The relative shift AG V + AF V is dependent upon u> via 
the atomic factors a n , a' n , and cr„; thus the lattice in- 
duced frequency shift may be modulated by tuning the 
lattice lasers to an appropriate wavelength, the so-called 
magic wavelength. With the intent of eliminating the 
Ml and E2 effects in lattice clocks which scale with the 
square-root of intensity (second term of Eq. (fT5"]) ). Ka- 
tori et al. Q recently defined the magic wavelength by 
the criteria Act = and proposed experimental deter- 
mination of it by analyzing, and minimizing to zero, the 
w-dependence of A£ v . We see here that such an exper- 
imental procedure actually absorbs the second term of 
Eq. as well; the definition of the magic wavelength 
is naturally extended to eliminate the net (vibrational 
and non-adiabatic) effects which scale with the square- 
root of intensity. For operation at this magic wavelength, 
the remaining terms which give rise to a frequency shift 
are the first terms of Eqs. (|14I15[) . 

AG + AF = 12h (-Aa! + IRAa) . 

An interesting result is that the non-adiabatic shift only 
affects blue-detuned (1=1) lattice clocks to within our 
approximations. 



We now estimate the size of the shifts AG and AF for 
the prospective Sr lattice clock operating at the 390 nm 
blue magic wavelength. We find a total shift 

AG + AF = Y2\a\I a (1.4 x 1(T 7 + 9.5 x 1CT 10 ) , (16) 

where oip = a$ = ct and the two terms on the right- 
hand-side are respective contributions. To obtain this 
result we have borrowed the ratio Aa'/|a| = —1.4 x 10~ 7 
from Ref. [|| and we have evaluated equations (|10I12[) 
using dipole matrix elements from Ref. [l5j and energies 
from Ref. [3 (a s = a P = -455 a.u., a s = 28000 a.u., 
up = 218000 a.u.). Assuming I x = I y = I z = Iq with 
a potential depth of 10 fiK (the corresponding potential 
depth-to-recoil energy ratio being A\a\Io/R = 14), we 
find the following fractional shifts to the clock frequency 

5v AG + AF lfi 18 

— = = 2.0 x 10~ 16 + 1.4 x 10" 18 . 

v AW 

We see that the non-adiabatic shift is on the order of 
10~ 18 , accounting for about 0.7% of the total shift. This 
is on the same order as the Ml contribution (0.8% Q), 
both of these effects being overshadowed by the E2 con- 
tribution. 

We note that for the assumed trap depth and lat- 
tice geometry, clock accuracy at the 10 -18 level can 
only be achieved if the combined parameter Aa'Io is 
known/controlled to better than 5%. This may prove 
difficult to attain experimentally. In Ref. two novel 
lattice geometries were conceived such as to effectively 
eliminate cither the E2 (case (II)) or the Ml (case (III)) 
frequency shift. For the blue Sr lattice clock, implemen- 
tation of the case (II) geometry may be necessary to real- 
ize accuracy at the 10 -18 level. In such a set-up, the non- 
adiabatic and Ml contributions are no longer masked by 
the E2 shift and are likely to constitute the leading un- 
certainty in the clock. 

It is also worth noting that the red-detuned 3-D lattice 
clock of Ref. [3 does not have the specific geometry 
considered above and thus may be susceptible to the non- 
adiabatic effects. 

If we treat the potential depth-to-recoil energy ratio as 
a fixed parameter, then the non-adiabatic shift is found 
to scale as R 2 . Thus we reason that the non-adiabatic 
effects may be enhanced for blue lattice clocks based on 
other atomic species. Considering the lattice clock can- 
didate Mg (see, e.g., Ref. [I3l); wc note that the atomic 
mass and energy spectrum necessitate a larger recoil en- 
ergy than the blue Sr lattice clock (i?M g ^ 7i?sr) and 
the resulting non-adiabatic shift could easily reach to the 
10 -16 level. We further note that the non-adiabatic ef- 
fects could potentially be enhanced for blue magic wave- 
lengths which are relatively close to atomic resonance, as 
is evident from the approximate relation <j n /a n ~ —1/S n , 
with S n being the energy detuning from resonance. 
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